HEAT EXCHANGER OPTIMIZATION PROCESS AND APPARATUS 



BACKGROUND OF THE INVENTION 

01 The energy required to transport natural gas through a pipeline is stored in the form of gas 
pressure. Therefore as natural gas is transported along a pipeline it undergoes a pressure loss. For 
any given pipeline and gas flow rate, the magnitude of the pressure loss is primarily dependent upon 
the pressure and temperature of the gas flowing through the pipeline. The lower the natural gas 
pressure and/or higher the temperature, the greater the pipeline pressure loss, and vice versa. 

02 In order to maximize natural gas throughput of a pipeline it is necessary to restore the 
pressure loss (i,e. energy consumed) that occurs as the natural gas travels down a pipeline, by 
compressing the gas at regular intervals along the pipeline. For any given compressor discharge gas 
pressure and gas flow rate, the power and thereby energy required to compress the gas is primarily 
dependent upon the compressor suction gas pressure and temperature (i.e. pipeline outlet gas 
pressure and temperature). The lower the suction gas pressure and/or higher the suction gas 
temperature, the greater the amount of power required to compress the natural gas, and vice versa. 

03 However natural gas undergoes an increase in temperature during the compression process. 
Depending on the length of intervening piping/pipeline and surrounding soil conditions, elevated 
compressor station discharge gas temperatures (i.e. pipeline inlet gas temperatures) result in lower 
suction gas pressures and/or higher suction gas temperatures at the downstream compressor station. 
For any given gas flow rate this results in an increase in the amount of power required to compress 
the natural gas at the downstream compressor and thereby results in an increase in the amount of 
energy required to transport the gas through the pipeline. 

04 Air cooled heat exchangers are utilized to lower the temperature of the natural gas to reduce 
compression power requirements and thereby reduce the amount of energy required to transport 
natural gas through a pipeline. 
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05 In the case of air cooled natural gas heat exchangers, heat transfer occurs between the warm 
gas and cool ambient air that is forced through the exchanger by powered fans. The amount by 
which the gas temperature is reduced is dependent upon heat exchanger design (i.e. total cooling 
surface area, and number, configuration, length and diameter of the cooling tubes, and overall heat 
transfer coefficient), heat exchanger inlet gas pressure and temperature, ambient air temperature, the 
specific heat capacity of the natural gas and air, and flow rate of the gas and air througji the 
exchanger. 

06 It is known in the art that the cooling of natural gas is sometimes not energy efficient and 
consequently, heat exchangers are often shut down and/or bypassed completely to avoid 
unnecessary energy costs. This invention provides an improvement upon prior art methods of 
optimizing energy savings during the operation of compressor stations and heat exchangers when 
transporting natural gas through pipelines. 

SUMMARY OF THE INVENTION 

07 The inventor has recognized that a pressure loss occurs as natural gas passes through a heat 
exchanger, which will offset the power and thereby energy savings afforded by cooling the gas. For 
any given inlet gas pressure and temperature, this pressure loss varies with the amount of gas flow 
through the heat exchanger. 

08 An increase in gas flow rate through a heat exchanger results in greater heat transfer (i.e. 
greater cooling of the total gas flow coming fi-om the upstream compressor), and vice versa. 
However an increase in gas flow rate also results in a greater pressure loss for the natural gas, and 
vice versa. 

09 Given the manner in which the amount of heat transfer (i.e. cooling) and heat exchanger 
pressure loss varies with the amount of gas flow through the exchanger, there is, for any given heat 
exchanger design (i.e. total cooling surface area, and number, configuration, length and diameter of 
the cooling tubes, overall heat transfer coefficient, and cooling fan power requirement ) and 
operating conditions (i.e. heat exchanger inlet gas pressure and temperature, ambient air 
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temperature, specific heat capacity of the natural gas and air, total flow rate of natural gas from 
upstream compressor(s), and air flow rate through the heat exchanger), an optimum gas flow rate 
through the exchanger that yields the maximum net compression power (and thereby energy) 
savings. 

10 There is therefore provided in accordance with an aspect of the invention, an electronic Heat 
Exchanger Management Device (HEMD) that determines the optimum balance between gas cooling 
and heat exchanger pressure loss, for any given operating condition. The HEMD then adjusts the 
gas flow rate through the exchanger accordingly, to yield the maximum net power (and thereby 
energy) savings obtainable from that exchanger. The HEMD can maintain the optimum balance 
between gas cooling and heat exchanger pressure loss and reduce the amount of energy required to 
transport a given volume of gas through a pipeline and thereby increase the transmission efficiency 
of the gas pipeline system. 
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1 1 According to a further aspect of the invention, there is provided a method of operating a heat 
exchanger on a natural gas transmission pipehne using a control algorithm. The control algorithm is 
the heart of the HEMD. The control algorithm is used to calculate the optimum position of the heat 
exchanger bypass valve to achieve the optimum gas flow rate through a exchanger based upon heat 
exchanger design (i,e. total cooling surface area, and number, configuration, length and diameter of 
the cooling tubes, overall heat transfer coefficient, and cooling fan power requirement) and 
operating conditions (i.e. heat exchanger inlet gas pressure and temperature, ambient air 
temperature, specific heat capacity of the natural gas and air, total flow rate of gas firom upstream 
compressor(s), and air flow rate through the heat exchanger). The operating conditions of the heat 
exchanger and upstream compressor are continually monitored for any changes by the HEMD and 
the control algorithm is utilized to determine, among other things, the optimum flow through the 
heat exchanger and optimxmi position of the heat exchanger bypass valve, and to initiate a command 
signal to the bypass valve operator to change the position of the bypass valve, accordingly. Upon 
receipt of a command signal the bypass valve will then close or open and thereby increase or 
decrease the gas flow rate through the exchanger, as required. The command signal to the bypass 
valve operator will continue until the optimum balance between gas coolmg and heat exchanger 
pressure loss (as calculated by the control algorithm) has been achieved and detected by the HEMD. 

12 The control algorithm is derived through computer modeling of the heat exchanger, 
associated upstream and downstream compressors, and the intervening pipeline system(s). The 
model is used to determine the impact of incremental changes of flow through the exchanger on 
pipeline compression power requirements, based on a number of varying operating parameters (i.e. 
total flow rate of natural gas fi-om the upstream compressor(s), air flow rate through the heat 
exchanger, heat exchanger inlet and outlet gas pressure and temperature, ambient air temperature, 
and power (and thereby energy) requirements of the upstream gas compressor and heat exchanger 
cooling fans). Where possible, the heat exchanger performance characteristics utilized in the model 
are based upon actual site measured parameters in order to improve the accuracy of the control 
algorithm. The results of this modehng form the basis for the control algorithm. 

13 The resulting control algorithm takes into account the physical characteristics of the actual 
upstream and downstream compressors (i.e. compressor power requirement, performance and 
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efficiency), the intervening pipeline systems (i.e. pipeline size, length, surrounding soil conditions, 
and heat conductivity) and heat exchanger design (i.e. total cooling surface area, and number, 
configuration, length and diameter of the cooling tubes, cooling air flow capacity, and overall heat 
transfer coefficient) and is therefore unique to each heat exchanger installation. The control 
algorithm also takes into account the current operating conditions of the heat exchanger (i.e. heat 
exchanger inlet and outlet gas pressure and temperature, ambient air temperature, the specific heat 
capacity of the natural gas and air, and air flow rate through the heat exchanger) and the upstream 
compressor(s) (i.e. total gas flow rate, discharge gas pressure and temperature). The control 
algorithm can among other things define the relationship between the temperature differential 
between the heat exchanger inlet gas and ambient air, and the pressure loss across the exchanger, 
that yields the optimum gas flow rate through the exchanger and thus the maximum net power (and 
thereby energy) savings. 

14 According to a further aspect of the invention, given the manner in which the control 
algorithm is derived, it can be utilized to compare current (i.e. real time) actual heat exchanger 
performance to predicted performance criteria (based on actual site measured baseline parameters). 
As a result the HEMD can among other things, be utilized to monitor the extent of internal or 
extemal exchanger fouling (i.e. heat exchanger performance) and to determine if and when the 
exchanger may require cleaning/de-fouling. 

BRIEF DESCRIPTION OF THE FIGURES 

15 There will now be described preferred embodiments of the invention, with references to the 
figures by way of illustration, in which like reference characters denote like elements, and in which: 

Fig. 1 is a schematic showing a natural gas pipeline with an upstream compressor station 
"A", heat exchanger and downstream compressor station "B"; 

Fig, 2A and 2B is a flow diagram illustrating method steps according to the invention; 

Fig. 3 is a schematic showing inputs and outputs of a first embodiment of a heat exchanger 
management device according to the invention; 

Fig. 4 is a schematic showing inputs and outputs of a second embodiment of a heat 
exchanger management device according to the invention; 
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Fig. 5 is a schematic showing inputs and outputs of a third embodiment of a heat exchanger 
management device according to the invention; 

Fig. 6 is a schematic shov^^ing inputs and outputs of a fourth embodiment of a heat exchanger 
management device according to the invention; 

Fig. 7 is a schematic showing inputs and outputs of a fifth embodiment of a heat exchanger 
management device according to the invention; 

Fig. 8 is a schematic showing inputs and outputs of a sixth embodiment of a heat exchanger 
management device according to the invention; 

Fig. 9 is a schematic showing inputs and outputs of a seventh embodiment of a heat 
exchanger management device according to the invention; 

Fig. 10 is a schematic showing inputs and outputs of a eighth embodiment of a heat 
exchanger management device according to the invention; 

Fig. 11 is a schematic showing inputs and outputs of a ninth embodiment of a heat 
exchanger management device according to the invention; 

Fig. 12 is a schematic showing inputs and outputs of a tenth embodiment of a heat 
exchanger management device according to the invention; 

Fig. 13 is a schematic showing inputs and outputs of a eleventh embodiment of a heat 
exchanger management device according to the invention; 

Fig. 14 is a graph showing net power savings through cooling by use of a method and 
apparatus according to the invention; 

Fig. 15 is a graph showing optimum amount of cooling obtained by operation of a method 
and apparatus according to the invention; 

Fig. 16 is a graph showing benefits of cooling at various ambient temperatures fi*om 
operation of a method and apparatus according to the invention; and 

Fig. 17 is a graph illustrating the basis for the control algorithm for optimizing power (and 
thereby energy savings), based on various temperature differentials between the heat exchanger 
inlet gas and ambient air, for a given gas flow rate. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS OF THE INVENTION 
16 In this patent document, the word comprising is used in its inclusive sense and does not 
exclude other elements being present. The indefinite article "a" before an element does not exclude 
other elements being present. 



17 The energy required to transport natural gas through a pipeline is stored in the form of gas 
pressure. Therefore as natural gas is transported along a pipeline it undergoes a pressure loss. The 
relationship of gas mass flow rate, pipeline inlet gas pressure and temperature to pipeline pressure 
loss can be defined as follows (based on the AGA Equation for pipeline pressure loss given fixlly 

turbulent pipeline flow); 
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18 Given that standard pressure and temperature conditions, and pipeline length, relative 
roughness and diameter and natural gas specific gravity are constants then the impact of changing 
gas mass flow rate, pipeline inlet gas pressure and temperature to pipeline pressure loss can be 
defined as follows; 
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19 A decrease in pipeline inlet natural gas pressure will result in an increase in pipeline 
pressure loss, and vice versa. Assuming that the gas mass flow rate, gas temperature and gas 
compressibility are essentially constant, it can be shown that lowering the gas pressure at the inlet of 
the pipeline by a given amount will result in an increase in pipeline pressure loss, as follows; 
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20 An increase in pipeline inlet natural gas temperature will increase pipeline pressure loss, and 
vice versa. Assuming that pipeline inlet gas pressure, gas mass flow rate and gas compressibility are 
essentially fixed, it can be shown that the amount of pipeline pressure loss is dependent upon the 
temperature of the gas flowing through the pipeline, as follows; 
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21 In order to maximize the natural gas throughput of a pipeline it is necessary to restore the 
pressure loss (i.e. energy consumed) that occurs as gas travels through a pipeline by compressing 
the gas at regular intervals along the pipeline. The power required to compress gas is a function of 
compressor gas mass flow rate, compressor efficiency, compressor suction gas pressure and 
temperature, discharge gas pressure, and gas composition/characteristics (i.e. adiabatic exponent, 
compressibility, and molecular weight), as follows; 
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22 A decrease in suction gas pressure and/or increase in discharge gas pressure will resuh in an 
increase in the amount of power required for compression, and vice versa. Assuming that natural 
gas mass flow rate, compressor efficiency, and gas composition/characteristics (i.e. adiabatic 
exponent, compressibility, and molecular weight) and suction gas temperature are essentially fixed, 
it can be shown that the amount of compression power required is dependent upon suction and 
discharge gas pressure, as follows; 

>^[(wrO/*a*«.] ^ 
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23 An increase in the compressor suction gas temperature will resuh in an increase in the 
amount of power required for compression, and vice versa. Assuming that natural gas mass flow 
rate, compressor efficiency, and gas composition/characteristics (i.e. adiabatic exponent, 
compressibility, and molecular weight) and suction and discharge gas pressure are essentially fixed, 
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it can be shown that the amount of compression power required is dependent upon suction gas 
temperature, as follows; 



Pvt/r = Pwr M * 
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24 Gas undergoes an increase in temperature during the compression process, as follows; 



T = r +AT 
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25 A decrease in suction gas pressure and/or increase in discharge gas pressure will result in a 
higher compressor discharge gas temperature and thereby a higher downstream pipeline inlet gas 
temperature, and vice versa. Assuming that compressor efficiency, natural gas 
composition/characteristics (i.e. adiabatic exponent) and suction gas temperature are essentially 
fixed, then it can be shown that compressor discharge gas temperature is dependent upon suction 
and discharge gas pressure, as follows: 
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26 An increase in suction gas temperature will resuU in a higher compressor discharge gas 
temperature and thereby a higher downstream pipeline inlet gas temperature, and vice versa. 
Assuming that compressor efficiency, natural gas composition/characteristics (i.e. adiabatic 
exponent) and suction and discharge gas pressure are essentially fixed, then it can be shown that 
compressor discharge gas temperature is dependent upon suction gas temperature, as follows; 
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27 Air cooled gas heat exchangers are utilized to lower the temperature of the compressor 
discharge gas as a means to reduce compression power and thereby energy requirements. 



28 hi the case of air cooled natural gas heat exchangers, heat transfer occurs between the warm 
natural gas and cool ambient air that is forced through the exchanger by powered fans. The amount 
by which the gas temperature is reduced is dependent upon heat exchanger design (i.e. total cooling 
surface area, and number, configuration, length and diameter of the cooling tubes, cooling air flow 
capacity, and overall heat transfer coefficient), heat exchanger inlet gas pressure and temperature, 
ambient air temperature, the specific heat capacity of the natural gas and air, mass flow rate of the 
gas and air through the exchanger. 

The heat transfer to the air equals; 

9 air ^^lirHE * ^ pair * ^airHEout " '^airamh {HEin)) 

which equals the heat transfer fix)m the gas; 
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However the actual amoimt of heat transfer that occurs between the natural gas and air is a function 
of the physical characteristics of the heat exchanger being utilized. Heat transfer by a heat 
exchanger can be defined as follows; 
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and 
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29 The amount of heat transfer that occurs between the natural gas and air, the temperature of the 
heat exchanger outlet gas, and the temperature of the heat exchanger outlet air, can be determined 
by solving the three equations simultaneously, as follows; 

HE* ^HEext * ^ * ^MTD = f&^asfjE * ^ pgas * tasHEin " ^gasHEout ) ~ ^tirHE * ^ pair * ^airHEout " ^airambiHEin) ) 

An increase in the natural gas mass flow rate through a heat exchanger will result in an increase in 
the overall heat transfer coefficient and thereby a greater heat transfer between the warm gas and 
cooling air, and vice versa. Given that the heat exchanger physical characteristics are constant and 
assuming that gas and air composition/characteristics (i.e. specific heat capacity, thermal 
conductivity and dynamic viscosity), heat exchanger inlet gas temperature, air mass flow rate, and 
ambient air temperature are essentially fixed, it can be shown that the overall heat transfer 
coefficient and amount of heat transfer are dependent upon the gas mass flow rate through the heat 
exchanger. 



Conversely, an increase in the air mass flow rate through a heat exchanger will result in an increase 
in the overall heat transfer coefficient and thereby a greater heat transfer between the warm natural 
gas and cooling air, and vice versa. Given that the heat exchanger physical characteristics are 
constant and assuming that gas and air composition/characteristics (i.e. specific heat capacity, 
thermal conductivity and dynamic viscosity), heat exchanger inlet gas temperature, gas mass flow 
rate, and ambient air temperature are essentially fixed, it can be shown that the overall heat transfer 
coefficient and amount of heat transfer are dependent upon the air mass flow rate through the heat 
exchanger. 
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30 The resultant temperature of the natural gas downstream of a heat exchanger, that results 
from the mixture of the heat exchanger outlet gas and any non-cooled bypass gas, can be 
determined as follows (Note bypass gas temperature equals heat exchanger inlet gas temperature); 

T gnsmixed ^ ([('^osconv ~ ^gasHB )* '^gasHEin J"^ [^gasHE * '^gasHEout j)'^ '^gascomp 

3 1 Gas passing through heat exchanger cooling tubes will undergo a pressure loss, as follows; 

^gasHB = 0M057 * ^"""^ X -^^"^^ * 



32 An increase in natural gas mass flow rate through a heat exchanger will result in a greater 
pressure loss across the heat exchanger, and vice versa. Assuming that gas composition (i.e. gas 
density) and heat exchanger tube friction factor, length, and diameter are essentially constant, it can 
be shown that the magnitude of the pressure loss is dependent upon the gas mass flow rate through 

the heat exchanger, as follows; 
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33 The effect lower natural gas temperature has on downstream compression conditions and 
thereby downstream compression power requirements depends on the length of the intervening 
piping/pipeline between the heat exchanger and downstream compressor. If the piping is short, any 
reduction in downstream compression power will result from a lower suction gas temperature (refer 
to Paragraphs 23 and 30). However the intervening pipeline may be of sufficient length such that 
the temperature of the gas leaving the pipeline is similar to the temperature of the gas had no 
upstream heat exchanger been present (i.e. equal to the groimd temperature). In this instance any 
reduction in downstream compression power will result from a lower pipeline pressure loss and 
thereby higher suction gas pressure that is the resuU of cooler gas entering the pipeline (refer to 
Paragraphs 20 and 22). For pipelines of an intermediate length any reduction in downstream 
compression power will result from a combination of lower suction gas temperature and higher 
suction gas pressure. 
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34 The majority of the effect lower natural gas temperature has on reducing downstream 
compression power requirements (refer to Paragraph 33) can be redirected to the upstream 
compressor. It can be shown that it is possible to reduce the discharge pressure and thereby power 
requirements of the upstream compressor such that after the gas has passed through the exchanger 
there is no net negative impact on the original downstream compression power requirements. The 
reduction in upstream compression power is possible whether or not the downstream compressor is 
located adjacent to or at some distance along a pipeline from the upstream compressor. 

35 The pressure loss associated with natural gas flowing through a heat exchanger will reduce 
the suction gas pressure at the downstream compressor and can result in an increase in compression 
power requirements (refer to Paragraphs 19, 22 and 32). As a result, depending on heat exchanger 
design (i.e. total cooling surface area, and number, configuration, length and diameter of the cooling 
tubes, and overall heat transfer coefficient), heat exchanger inlet gas pressure and temperature, and 
ambient air temperature, the specific heat capacity of the natural gas and air, and mass flow rate of 
the gas and air through the exchanger, the coolmg of the natural gas may in some instances result in 
an increase in the downstream compression power. This can occur when the pressure loss which 
results from the gas flowing through the heat exchanger, and the consequent energy loss, is greater 
than the energy savings afforded from the lower gas temperature 

36 It can be shown that the rate of decrease m compression power savings resulting from a 
reduction in gas cooling, that results from decreasing the gas mass flow rate through a heat 
exchanger from 100 percent flow, is initially much less than the rate of decrease in compression 
power losses resulting from the consequent reduction in pressure loss (refer to Paragraphs 19, 20, 
22, 23, 29, 30 and 32). 

37 It can be shown that it is possible to calculate the optimum natural gas mass flow rate which 
will yield the maximum net compression power savings. In simple terms the optimum gas mass 
flow rate will occur at the point at which the rate of decrease in compression power savings (as the 
gas mass flow rate through the heat exchanger is reduced and thereby the amount of gas cooling is 
reduced) equals the rate of decrease in compression power losses (resulting from the consequent 
reduction in pressure loss across the heat exchanger). It is at this point that the greatest difference 
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exists between compression power savings and losses. Further gas mass flow rate reduction will 
diminish the net savings because the rate of decrease in compression power savings becomes 
increasingly greater than the rate of decrease in compression power losses beyond this point. 

38 The optimum natural gas mass flow rate will vary for every heat exchanger installation. In 
particular the optimum gas mass flow rate through a heat exchanger will vary depending on heat 
exchanger design (i.e. total cooling surface area, and number, configuration, length and diameter of 
the cooling tubes, and overall heat transfer coefficient), heat exchanger inlet gas pressure and 
temperature, ambient air temperature, the specific heat capacity of the natural gas and air, total gas 
mass flow rate from upstream compressor, air mass flow rate through the heat exchanger. 

39 This invention provides an improvement upon prior art methods of optimizing energy costs 
by having the ability to determine the optimum natural gas mass flow rate through heat exchangers 
and control gas mass flow rate (when utilized in conjunction with an automated heat exchanger 
bypass valve), and thereby maintain the optimum amount of gas mass flow rate through heat 
exchangers, under varying operating conditions. This invention can continually maximize net 
compression power (and thereby net energy) savings afforded by heat exchangers on a real time 
basis, by continuously managing the gas mass flow rate through heat exchangers, based on variables 
such as the total gas mass flow rate from the upstream compressor, air mass flow rate through the 
heat exchanger, heat exchanger inlet gas pressure and temperature, ambient air temperature, and 
fouling resistance of the heat exchanger. 

40 Referring to Fig. 1, there is shown a natural gas pipeline 10 upon which is located an 
upstream compressor 12. After the upstream compressor 12, the gas pipeline 10 continues through a 
heat exchanger bypass valve 14, through a station gas flow meter 16 and on to a downstream 
compressor 18. The bypass valve 14 will typically include a conventional equalization bridle piping 
assembly and bypass valve operator (not shown). Closing of the bypass valve 14 diverts gas flow in 
the gas pipeline 10 to the heat exchanger piping 20 which passes through a heat exchanger gas flow 
meter 24 and heat exchanger 22 and then retums to the gas pipeline 10 downstream of the bypass 
valve 14. All of the elements 10 to 24 are conventional. 
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41 Various sensors/transmitters are provided on these elements, depending on the embodiment of 
the invention that is used. The sensors make observations of operating parameters that are 
characteristic of gas flow through the compressor station "A" and transmit the observation to the 
heat exchanger management device (HEMD) hereinafter referred to as the "controller 46". 
Hereinafter, these devices will be referred to as transmitters. A temperature transmitter 36 (output is 
TgasHEin) IS providcd on the inlet piping of the heat exchanger 22 to measure the heat exchanger inlet 
gas temperature. A temperature transmitter 38 (output is. TgasHEom) is provided on the outlet piping 
of the heat exchanger 22 to measure the heat exchanger outlet gas temperature. Temperature 
transmitters 29 (output is TairHEout) and 28 (output is Tairamb(HEin)) and pressure transmitter 33 (Paimmb) 
are provided near the heat exchanger 22 to measure the heat exchanger outlet air temperature and 
ambient air temperature and pressure, respectively. A temperature transmitter 44 (output is Tgasmixed) 
is provided downstream of the heat exchanger piping 20 and gas pipeline 10 junction to measure the 
temperature of the recombined gas flow. A pressure transmitter 34 (output is PgasHEin) is provided on 
the inlet piping of the heat exchanger 22 to measure the inlet gas pressure. A differential pressure 
transmitter 30 (output is APgasHE) is provided on either the bridle piping assembly of the bypass 
valve 14 or between the heat exchanger inlet and outlet headers (i.e. on either side of the heat 
exchanger 22) to measure the pressure loss across the heat exchanger 22. A heat exchanger gas 
flow meter 24 complete with differential pressure transmitter 32 (output is APgasmeter(HE)) in 
conjunction with pressure transmitter 34 (output is PgasHEin) and temperature transmitter 36 (output 
is TgasHEin) IS provided on the exchanger piping 20 to measure the heat exchanger gas flow rate. A 
station gas flow meter 16 complete with pressure transmitter 40 (output is Pgasmeter(station)) and 
differential pressure transmitter 42 (output is APgasmeter(station)) in conjunction with temperature 
transmitter 44 (output is Tgasmixed) is provided on the gas pipeline 10 to measure the total compressor 
station gas flow rate. A valve position transmitter 26 is provided on the heat exchanger bypass valve 
14 operator to provide bypass valve position information (output is BVposUion). Outputs from the 
transmitters 26 to 44 provided to controller 46, are described in more detail in relation to Figs. 3 to 
13. The controller 46 is preferably located adjacent to the heat exchanger bypass valve operator, 
where a suitable 12 to 24 VDC power supply is usually available. The controller 46 provides an 
output command signal (e.g. 12 to 24 VDC) to the heat exchanger bypass valve operator to either 
open (i.e. BVopen) or close (i.e. BVdosed) the bypass valve 14, as required. Elements 10, 12, 14, 16, 
20, 22, 24 and 26 to 46 together constitute a compressor station "A". A further compressor station 
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"B" downstream of compressor station "A" includes the gas pipeline 10 and downstream 
compressor 18. 

42 In a method for controlling natural gas cooling in a gas pipeline, the flow of gas through a heat 
exchanger at a compressor station, preferably including the effects of components upstream and 
downstream from the compressor station, is first modeled. The results of the modeling are then 
analyzed based upon the principles described previously in Paragraphs 17 to 38. Based on the 
application of these principles and current operating conditions of the compressor station a balance 
is found between gas cooling and heat exchanger pressure loss that results in an improvement of 
energy savings, preferably close to an optimal level. Finally, a bypass valve on the gas pipeline is 
operated to divert an amount of gas into the exchanger that achieves the desired balance. 

43 Applying the principles described previously in Paragraphs 17 to 38 and balancing the 
natural gas cooling and heat exchanger pressure loss based on the application of these principles and 
current operating conditions of the compressor station is carried out by a control algorithm that is 
preferably embodied within the controller 46, which may be a general purpose flow computer 
programmed according to this patent description. The following steps illustrate the development of 
the model and control algorithm with reference to Fig. 2A and 2B. 

44 Utilizing a flow analysis software (i.e. Simulation Science PRO/II, HYSYS or similar), 
create a model of compressor and heat exchanger and upstream and/or downstream pipeline(s), 
compressor(s), and heat exchanger(s) (step 50). 

45 Run the model based on the current actual operating conditions (i.e. gas flow rate, gas 
temperature and pressure conditions) of heat exchanger to be optimized. Determine whether or not 
the calculated temperature and pressure conditions upstream and downstream of heat exchanger 
reflect the actual conditions within an acceptable level of accuracy? If not, then modify the model as 
required (step 52). If yes, then proceed to step 54. 

46 Determine the actual temperature differential between the heat exchanger inlet gas and 
ambient air for the average Spring, Summer, Fall, Winter, and peak high and peak low ambient air 
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temperature day, based on operating records for the heat exchanger being modeled. Record this 
information (step 54). 

47 Determine the actual average gas flow rate and peak high and low gas flow rates through the 
compressor and the corresponding compressor discharge gas pressure at these flow rates, based on 
operating records for the compressor associated with heat exchanger being modeled. Record this 
information (step 56). 

48 Determine the percentage of compressor gas flow typically diverted through the heat 
exchanger and the corresponding pressure loss across the exchanger, given the actual average gas 
flow rate and peak high and peak low gas flow rates through the compressor, based on operating 
records for the heat exchanger being modeled. Record this information (step 58). 

49 Determine the percentage of the total number of heat exchanger cooling fans that typically 
would be utilized to cool the gas during the average Spring, Summer, Fall, Winter and peak high 
and peak low ambient air temperature day, based on operating records for the heat exchanger being 
modeled. Record this information (step 60). 

50 Run the model using the compressor gas flow rate, compressor discharge gas pressure and 
the temperature differential between the heat exchanger inlet gas and ambient air, corresponding to 
the average Spring day, with 100 % of the compressor gas flow rate diverted through the exchanger 
(i.e. 0 % exchanger bypass flow) and all of the cooling fans operating. Record the following 
information from the output data (step 62); 

percentage of gas flow through the heat exchanger 
pressure loss across the heat exchanger 

compressor suction and discharge gas pressure and temperature for all of the compressors in 
the model 

power requirements for all of the compressors in the model 

incremental change in compression power from the baseline power requirement (identified 
in the operating records) for all of the compressors in the model 
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51 Vary the gas flow rate through the heat exchanger from 100 % to 0 % of the total 
compressor gas flow rate, in 10 % increments, recording the same data as stated above in step 62, 
for each gas flow rate increment (step 64). 

52 Plot the heat exchanger gas flow rate percentage along the "X" axis and the total incremental 
change in compression power (from all of the compressors in the model) along the "Y" axis (step 

66). 

53 From the graph plotted in step 66 determine the minimum power required to transport the 
given gas flow rate through the pipeline(s). This point coincides with the optimum heat exchanger 
gas flow rate (step 68). 

54 Determine the optimum heat exchanger gas pressure differential associated with the 
optimum exchanger gas flow rate and record this information together with the corresponding 
temperature differential between the heat exchanger inlet gas and ambient air (step 70). 

55 Repeat steps 62 through 70 for various temperature differentials between the heat exchanger 
inlet gas and ambient air to obtain at least five other optimum operating points (i.e. repeat for the 
average Summer, Fall, Winter and peak high and peak low ambient air temperature days) (step 72). 

56 Plot the temperature differential between the heat exchanger inlet gas and ambient air along 
the "X" axis and the optimum heat exchanger gas pressure differential along the "Y" axis for the six 
optimum operating points produced after completing step 72 (step 74). 

57 Determine the mathematical relationship between the temperature differential (between the 
heat exchanger inlet gas and ambient air) and the optimum heat exchanger gas pressure differential 
based on the graph constructed in step 74. The mathematical equation for the line joining the 
different data points results in a basic heat exchanger control algorithm for the modeled exchanger 
(step 76). 
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58 Repeat steps 62 to 76 for at least two different gas flow rates (i.e. peak high and peak low 
gas flow rates through the compressor) to determine the sensitivity of the heat exchanger control 
algorithm to changes in flow rates. This information can then be utilized to modify the control 
algorithm to accoimt for changes in gas flow rate (step 78). 

59 Determine whether or not the temperature and pressure conditions upstream and 
downstream of heat exchanger predicted upon implementation of the control algorithm reflect the 
actual conditions within an acceptable level of accuracy? If not, then modify the control algorithm 
as required (step 80). If yes, then no further steps are required. 

60 When compared with traditional heat exchanger operation at a given gas flow rate, 
optimizing heat exchanger operation results in either the transportation of more gas for the same 
compression power (i.e. energy) input or the transportation of the same amount of gas for less 
power (i.e. energy) input. 

61 Heat exchanger optimization can be utilized to either decrease the amount of compression 
power required from the compressor 12 immediately upstream of the exchanger 22, or it can be 
utilized to lower the compression power requirements of the downstream compressors 18 (i.e. 
compressor 18). In the latter instance heat exchanger optimization can lower the amount of 
compression power required for the next two to three downstream compressors and thereby can 
result in greater energy savings. Whether the heat exchanger optimization will reduce the upstream 
and/or downstream compression power requirements depends upon how the pipeline system is 
modeled. In order to create a control algorithm that lowers the compression power requirements of 
the downstream pipeline system, the compressor immediately upstream of the exchanger should be 
modeled to control on constant compressor discharge gas pressure. In order to create a control 
algorithm that decreases the amount of compression power required for the compressor inamediately 
upstream of the exchanger, the upstream compressor station should be modeled to control on 
constant station discharge gas pressure. 

62 In order to determine the impact of heat exchanger optimization it is necessary to first 
establish a baseline. This can be done by reviewing the operating history for the upstream and/or 
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downstream compressor(s) and comparing the power requirements versus gas throughput ratio 
before heat exchanger operation has been optimized, to the same ratio after optimization has 
occurred (note the impact of any re-cahbration of instrumentation involved with these 
measurements has to be assessed prior to derivmg these ratios). This requires monitoring the gas 
flow rate through the upstream and/or downstream compressor(s) and the fuel (i.e. energy) 
consumption of the compressor drivers (e.g. gas turbines). 

63 The inputs to and outputs from a controller 46 according to the invention (see Fig. 1) are 
shown in Figs. 3 to 13. In all ten embodiments shown in Figs. 3 to 13, the controller 46 controls gas 
flow through a heat exchanger. Transmitters 30, 34, 36, 38, 28, 29 (optional), and 33 (optional) are 
provided with all embodiments. The APgasHE , PgasHEin 5 TgasHEin 5 TgasHEout ? Tairamb(HEin) ? TairHEout 
(optional), and Pairamb (optional) signals from these transmitters are acquired and routed into the 
controller 46. The controller 46 monitors the actual heat exchanger gas pressure differential 
(APgasHE) and compares this to the optimum heat exchanger gas pressure differential based upon a 
pre-programmed control algorithm. If the actual heat exchanger gas pressure differential (APgasHE) 
varies from the optimum heat exchanger gas pressure differential (as calculated by the control 
algorithm) the controller 46 sends a command signal (e.g. 12 to 24 VDC) to the operator of the heat 
exchanger bypass valve 14 to either open (BVopen) or close (BVdose) the bypass valve 14, as 
required. The command signal to move the bypass valve 14 will continue until the gas pressure 
differential (APgasHs) equals the optimum differential pressure calculated by the control algorithm. 
The controller 46 is based on commercially available flow computer technology together with 
added circuitry to accommodate the control algorithm which is calibrated/programmed with actual 
heat exchanger operating conditions and gas flow rate at the time of commissioning (Data in). The 
accuracy with which the control algorithm optimizes heat exchanger operation depends upon how 
closely the algorithm reflects the actual operating conditions of the heat exchanger. However, heat 
exchanger operating conditions can change with time. In order to ensure the accuracy of the control 
algorithm, programming updates for the algorithm are periodically required. In order to provide 
these programming updates it is first necessary to retrieve information regarding controller 46 status 
and actual heat exchanger operating parameters. These programming uploads and information 
downloads are possible through interfacing with either a local or remote upload/download devices 
(e.g. computer) via the controller 46 data input (Data in) and data output (Data out)? respectively. 
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64 In a first embodiment shown in Fig. 3, transmitters are provided as discussed in Paragraph 
60. Heat exchanger gas flow measurement (QgasstdHEest) and periodic compressor station gas flow 
measurement (Qgasstdstationest), are estimated in the manner described in Paragraphs 76B and 76M, 
respectively. The accuracy of these estimated gas flow measurements is contingent upon there being 
no change in internal fouling of the heat exchanger(Rfouiing). This is because a change in internal 
fouling (Rfouiing) would affect the heat exchanger gas pressure differential (APgasHE) in a similar 
fashion as a change in gas flow rate (i.e. an increase in internal fouUng would result in an increase 
the pressure differential as would an increase in gas flow rate), hi this embodiment both the 
estimated heat exchanger gas flow measurement (QgasstdHEest) and compressor station gas flow 
measurement (Qgasstdstationest) calculations utihze heat exchanger gas pressure differential (APgasHE) 
data. Should a change in internal foulmg (Rfouiing) occur then these esthnated gas flow 
measurements would be compromised. 

65 hi a second embodiment shown in Fig. 4, transmitter 26 is provided together with the 
transmitters included in embodunent one. The BVposition signal from this transmitter is acquked and 
routed into the controller 46. The second embodiment includes all of the functionality of 
embodiment one. Actual heat exchanger bypass valve gas flow measurement (Qgasstdbypassactuai) is 
calculated in the manner described in either Paragraph 76H or 761. Estimated compressor station 
gas flow measurement (Qgasstdstationest) is calculated based on estimated heat exchanger gas flow 
measurement (QgasstdHEest) and actual heat exchanger bypass valve gas flow measurement 
(Qgasstdbypassactuai), in the manner described in Paragraph 76N. The accuracy of the estimated gas flow 
measurements is contingent upon there being no change in internal fouling of the heat exchanger 
(Rfouiing), as with embodiment one and would be compromised should a change in internal fouling 
(Rfouiing) occur. It is possible to detect whether or not internal and/or external fouling of the heat 
exchanger (Rfouiing) has occurred if actual compressor station gas flow measurement (Qgasstdstationactuai) 
data is available (i.e. from outside sources m the case of embodiment two) over the same time 
period as gas flow measurement data recorded by the controller 46. The actual heat exchanger gas 
flow measurement (QgasstdHEactuai) can be calculated based on actual heat exchanger bypass valve gas 

flow measurement (Qgasstdbypassactuai) a^d compressor station gas flow measurement (Qgasstdstationactuai) 

in a manner described in Paragraph 76D. Any change in internal and/or external fouling of the heat 
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exchanger (RfouUng) can be detected and quantified based on original and actual heat exchanger gas 
flow measurement (QgasstdHEong and QgasstdHEactuai, respectively), in the manner described in 
Paragraphs 76T, 76U, 76V, 76W, 76X, 76Y, 76Z, 76AA, 76BB, and 76CC. This information can 
then be utilized to modify the control algorithm, if necessary. 

66 In a third embodiment shown in Fig. 5, transmitter 32 is provided together with the 
transmitters included in embodiments one and two. The heat exchanger gas flow meter pressure 
differential (APgasmeter(HE)) signal from this transmitter is acquired and routed into the controller 46. 
The third embodiment includes all of the functionality of embodiments one and two. Actual heat 
exchanger gas flow measurement (QgasstdHEactuai) is calculated in the manner described in Paragraph 
76C, Actual compressor station gas flow measurement (Qgasstdstationactuai) is calculated based on actual 
heat exchanger gas flow measurement (QgasstdHEactuai) and heat exchanger bypass valve gas flow 
measurement (Qgasstdbypassactuai), in the manner described in Paragraph 76Q. The accuracy of the 
actual gas flow measurement is not affected by a change in internal fouling of the heat 
exchanger(Rfouimg). Any change in internal and/or external fouling of the heat exchanger (RfouUng) 
can be detected and quantified based on original and actual heat exchanger gas flow measurement 
(QgasstdHEong and QgasstdHEactuah respectively), as described in embodiment two. As with embodiment 
two, this information can then be utihzed to modify the control algorithm, if necessary. 

67 In a fourth embodiment shown in Fig. 6, transmitters 42 and 40 are provided together with 
the transmitters included in embodiments one and two. The compressor station gas flow meter 
pressure differential (APgasmeter(station)) and compressor station gas flow meter pressure (Pgasmeter(station)) 
signals from this transmitter are acquired and routed into the controller 46. The fourth embodiment 
includes all of the functionaUty of embodiments one and two. Actual compressor station gas flow 
measurement (Qgasstdstationactuai) IS Calculated in the manner described in Paragraph 76P. Actual heat 
exchanger gas flow measurement (QgasstdHEactuai) is calculated based on actual compressor station gas 
flow measurement (Qgasstdstationactuai) and heat exchanger bypass valve gas flow measurement 
(Qgasstdbypassactuai), in the maimer described in Paragraph 76D. As with embodiment three, the 
accuracy of the actual gas flow measxirement is not affected by a change in internal fouling of the 
heat exchanger(Rfouiing). Any change in internal and/or external fouling of the heat exchanger 
(RfouUng) can be detected and quantified based on original and actual heat exchanger gas flow 
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measurement (QgasstdHEorig and QgasstdHEactuai, respectively), as described in embodiment two. As with 
embodiment two, this information can then be utiUzed to modify the control algorithm, if necessary. 

68 In a fifth embodiment shown in Fig. 7, transmitter 44 is provided together with the 
transmitters included in embodiment one. The gas temperature downstream of the heat exchanger 
piping and heat exchanger bypass piping junction (Tgasmixed) signal from this transmitter is acquired 
and routed into the controller 46. The fifth embodiment includes all of the fimctionality of 
embodiment one. Estimated heat exchanger gas flow measurement (QgasstdHEactuai) is calculated in 
the manner described in Paragraph 76B, as with embodiment one. Estimated heat exchanger bypass 
valve gas flow measurement (Qgasstdbypassest) is calculated based on estimated heat exchanger gas 
flow measurement (QgasstdHEest), in the manner described in Paragraph 76G. Estimated compressor 
station gas flow measurement (Qgasstdstationest) is calculated based on estimated heat exchanger gas 
flow measurement (QgasstdHEest) and gas temperature downstream of heat exchanger piping and heat 
exchanger bypass piping junction (Tgasmixed), in the manner described in Paragraph 760. The 
accuracy of the estimated gas flow measurements is contingent upon there being no change in 
internal fouling of the heat exchanger (Rfouiing), as with embodiment one and would be compromised 
should a change in internal fouling (Rfouiing) occur. It is possible to detect whether or not internal 
and/or external fouling of the heat exchanger (RfouUng) has occurred if actual compressor station gas 
flow measurement (Qgasstdstationactuai) data is available (i.e. from outside sources in the case of 
embodiment five) over the same time period as the gas flow measurement data recorded by the 
controller 46. The actual heat exchanger gas flow measurement (QgasstdHEactuai) can be calculated 
based on actual compressor station gas flow measurement (Qgasstdstationactuai) and estimated 
compressor station gas flow measurement (Qgasstdstationest) in a manner described in Paragraph 76GG. 
Any change in internal and/or external fouling of the heat exchanger (Rfouiing) can be detected and 
quantified based on original and actual heat exchanger gas flow measurement (QgasstdHEorig and 
QgasstdHEactuai^ respectively), in the manner described in Paragraphs 76T, 76U, 76V, 76W, 76X, 76Y, 
76Z, 76AA, 76BB, and 76CC. This information can then be utilized to modify the control 
algorithm, if necessary. 

69 In a sixth embodiment shown in Fig. 8, transmitter 32 is provided together with the 
transmitters included in embodiments one and five. The heat exchanger gas flow meter pressure 
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differential (APgasmeter(HE)) signal from this transmitter is acquired and routed into the controller 46. 
The sixth embodiment includes all of the functionality of embodiments one and five. Actual heat 
exchanger gas flow measurement (QgasstdHEactuai) is calculated in the manner described in Paragraph 
76C. Actual heat exchanger bypass valve gas flow measurement (Qgasstdbypassactuai) is calculated based 
on actual heat exchanger gas flow measurement (QgasstdHEactuai) and gas temperature downstream of 
heat exchanger piping and heat exchanger bypass piping junction (Tgasmixed), in the manner 
described in Paragraph 76K. Actual compressor station gas flow measurement (Qgasstdstationactuai) is 
calculated based on actual heat exchanger gas flow measurement (QgasstdHEactuai) and gas temperature 
downstream of heat exchanger piping and heat exchanger bypass piping junction (Tgasmixed), in the 
manner described in Paragraph 76R. The accuracy of the actual gas flow measurement is not 
affected by a change in internal fouling of the heat exchanger(Rfouiing)* Any change in internal 
and/or external fouling of the heat exchanger (Rfouiing) can be detected and quantified based on 
original and actual heat exchanger gas flow measurement (QgasstdHEong and QgasstdHEactuab 
respectively), as described in embodiment five. As with embodiment five, this information can then 
be utilized to modify the control algorithm, if necessary. 

70 hi a seventh embodiment shown in Fig. 9, transmitters 42 and 40 are provided together with 
the transmitters included in embodiments one and five. The compressor station gas flow meter 
pressure differential (APgasmeter(station)) and compressor station gas flow meter pressure (PgasmeterCstation)) 
signals from this transmitter are acquired and routed into the controller 46. The seventh embodiment 
includes all of the functionality of embodiments one and five. Actual compressor station gas flow 
measurement (Qgasstdstationactuai) is Calculated in the manner described in Paragraph 76P. Actual heat 
exchanger bypass valve gas flow measurement (Qgasstdbypassactuai) is calculated based on actual 
compressor station gas flow measurement (Qgasstdstationactuai) and gas temperature downstream of heat 
exchanger piping and heat exchanger bypass piping junction (Tgasmixed)^ in the maimer described in 
Paragraph 76L. Actual heat exchanger gas flow measurement (QgasstdHEactuai) is calculated based on 
actual compressor station gas flow measurement (Qgasstdstationactuai) and gas temperature downstream 
of heat exchanger piping and heat exchanger bypass piping junction (Tgasmixed), in the manner 
described in Paragraph 76F. As with embodiment six, the acciu-acy of the actual gas flow 
measurement is not affected by a change in internal fouling of the heat exchanger(Rfouiing). Any 
change in internal and/or external fouling of the heat exchanger (Rfouimg) can be detected and 
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quantified based on original and actual heat exchanger gas flow measurement (QgasstdHEong and 
QgasstdHEactuab respectively), as described in embodiment five. As with embodiment five, this 
information can then be utilized to modify the control algorithm, if necessary. 

71 In an eighth embodiment shown in Fig. 10, transmitter 26 is provided together with the 
transmitters included in embodiment five. The heat exchanger bypass valve position (BVposition) 
signal firom this transmitter is acquired and routed into the controller 46. The eighth embodiment 
includes all of the fiinctionality of embodiments one, two, and five. As with embodiment two, actual 
heat exchanger bypass valve gas flow measurement (Qgasstdbypassactuai) is calculated in the manner 
described in either Paragraph 76H or 761. Actual compressor station gas flow measurement 
(QgasstdstationactuaO is Calculated based on actual heat exchanger bypass valve gas flow measurement 
(Qgasstdbypassactuai) and gas temperature downstream of heat exchanger piping and heat exchanger 
bypass piping junction (Tgasmixed), in the manner described in Paragraph 76S. The accuracy of the 
actual gas flow measurement is not affected by a change in internal fouling of the heat 
exchanger(Rfouiing). Any change in internal and/or external fouling of the heat exchanger (Rfouiing) 
can be detected and quantified based on original and actual heat exchanger gas flow measurement 
(QgasstdHEong and QgasstdHEactual* respectively), as described in embodiments two and five. As with 
embodiments two and five, this information can then be utilized to modify the control algorithm, if 
necessary. 

72 In an ninth embodiment shown in Fig. 11, transmitter 32 is provided together with all of the 
transmitters included in embodiment eight. The heat exchanger gas flow meter pressure differential 
(APgasmeter(HE)) sigual from this transmitter is acquired and routed into the controller 46. The ninth 
embodiment includes all of the fiinctionality of embodiments one, two, three, five, six and eight. As 
with embodiment three, actual heat exchanger gas flow measurement (QgasstdHEactuai) is calculated in 
the manner described in Paragraph 76C. This embodiment provides two independent methods for 

calculating QgasstdHEactuai , Qgasstdbypassactuai j Qgasstdstationactual and thereby APgasmeter(HE) > UhEj and Rfouiing 

, which allows for cross checking of these calculations. This provides a greater degree of accuracy 
when modifying and reprogramming the control algorithm to obtain the optimum balance between 
gas cooling and heat exchanger pressure loss. 
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73 In an tenth embodiment shown in Fig. 12, transmitters 42 and 40 are provided together with 
all of the transmitters included in embodiment eight. The compressor station gas flow meter 
pressure differential (APgasineter(station)) £uid compressor station gas flow meter pressure (Pgasmeter(station)) 
signals from this transmitter are acquired and routed into the controller 46. The tenth embodiment 
includes all of the functionality of embodiments one, two, four, five, seven and eight. As with 
embodiment four, actual compressor station gas flow measurement (Qgasstdstationactua]) is calculated in 
the manner described in Paragraph 76P. This embodiment provides two independent methods for 

calculating QgasstdHEactual , Q gasstdbypassactual ? Qgasstdstationactua! and thereby APgasmeter(HE) , Uhe, and Rfouling 

, which allows for cross checking of these calculations. This provides a greater degree of accuracy 
when modifying and reprogramming the control algorithm to obtain the optimum balance between 
gas cooling and heat exchanger pressure loss. 

74 hi an eleventh embodiment shown in Fig. 13, transmitters include all of those found in 
embodiment eight and nine. The eleventh embodiment includes all of the functionality of 
embodiments one, two, three, four, five, six, seven, eight, nine and ten. This embodiment provides 
three independent methods for calculating QgasstdHEactual , Qgasstdbypassactuai , Qgasstdstationactuai and thereby 
APgasmeter(HE) , Uhe, and Rfouling , which allows for multiple cross checking of these calculations. 
This provides a greater degree of accuracy when modifying and reprogramming the control 
algorithm to obtain the optimum balance between gas cooUng and heat exchanger pressure loss. 

75 The following inputs and outputs are utiUzed and/or generated by the control algorithm; 

A) TgasHEin " heat exchanger inlet natural gas temperature, 

B) TgasHEout ~ heat exchanger outlet natural gas temperature, 

C) Tairamb(HEin) ~ ambient air temperature (i.e. heat exchanger inlet air temp),®K 

D) Pairamb - ambient air pressure, kPa(a) (optional) 

E) Tgasmixed " natural gas temperature resulting from the mixture of the heat exchanger outlet 
gas (TgasHEout) and the non-cooled bypass gas (TgasHEin), ""K 

F) PgasHEin - heat exchanger inlet natural gas pressure, kPa(a) 

G) APgasHE ~ natural gas pressure differential across heat exchanger, kPa(g) 
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H) APgasmeter(HE) " natural gas pressure differential across heat exchanger gas flow meter, kPa(g) 

I) APgasmeter(station) " natural gas pressure differential across compressor station gas flow meter, 
kPa(g) 

J) Pgasmeter(station) - natural gas prcssuTe at compressor station gas flow meter, kPa(a) 

K) TairHEout ~ heat exchanger outlet au: temperature,®K (optional) 

L) B Vposition -heat exchanger bypass valve actual position 

M) BVopen - heat exchanger bypass valve open command signal 

N) BVciosed - heat exchanger bypass valve closed command signal 

O) Power Supply - power supply to HEMD (12 to 24 VDC) 

P) Data out - HEMD operating data (historical and/or real time) for download to local and/or 
remote monitoring device(s). Depending on the type of sensors/transmitters available to the HEMD, 
data output can include some or all of the following information; TgasHEin , TgasHEom , Tairamb(HEm) » 

Pairamb , Tgastirixed , PgasHEin , APgasHE , APgasmeter(HE) Or APgasmeter(station) (or Calculated gaS floW rate 
"QgasstdHE" or "Qgasstdstation")* Pgasmeter(station) , TairHEout , BVposition (or Calculated AfiVopen), Calculated 

internal and external fouling heat transfer coefficients, power supply voltage, current control 

algorithm, and the time and date the data is recorded. 

Q) Data in - Calibration/programming data from local and/or remote upload device(s) to re- 
calibrate and update the HEMD software. Depending on the type of sensors/transmitters available to 
the HEMD, data input can include some or all of the following parameters; TgasHEin , TgasHEout , 

Tairamb(HEin) » Pairamb » Tgasmixed > PgasHEin 9 APgasHE > APgasmeter(HE) OV APgasmeter(station) 9 Pgasmeter(station) > 

TairHEout , BVposition , currcut internal and external fouling heat transfer coefficients, power supply 
voltage, updated control algorithm, and time and date the data is entered. 

76 The control algorithm may use some or all of the following calculations: 

A) Given ongoing PgasHEin » and TgasHEin data is available to the HEMD, then the natural gas 
specific gravity (SGgas) and compressibility factor (Zgas) are calculated based upon the on going data 
and the natural gas properties using the Benedict Webb Rubin (i.e. BWR) equation of state. 

B) Given ongoing APgasHE > PgasHEin , TgasHEin , Zgas 9 and SGgas data is available to the HEMD and 

that APgasHEorig 9 PgasHEinorig , TgasHEinorig 9 Zgasorig , SGgasorig > and QgasstdHEorig haS been previoUSly 

recorded when 100% of the natural gas flow was diverted through the heat exchanger during the 
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initial calibrationyprogramming of the HEMD, then the estimated gas flow rate through the heat 
exchanger (QgasstdHEest) can be calculated as follows; (Note the accuracy of this calculation is 
dependent upon the heat exchanger fouling resistance remaining constant) 
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C) Given on going APgasmeter(HE) , PgasHEin , TgasHEin , Zgas , SGgas , and heat exchanger gas flow 
orifice meter opening area (i.e. AorificeHE), data is available to the HEMD (Note Yonfice and Corifice are 
based on the relationship of AorificeHE to the cross-sectional area of the piping associated with the 
orifice meter and can be obtained from lookup tables), then the actual natural gas flow rate through 
the heat exchanger (QgasstdHEactuai) can be calculated as follows; 
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D) Alternatively to Calculations "C", "E" and "F", given on going Qgasstdbypassactuai (from 
Calculation "H" or "I", as appropriate or "L") and Qgasstdstationactuai (from Calculation "P", or "S"), 
then the actual natural gas flow rate through the heat exchanger (QgasstdHEactuai) can be calculated as 
follows; 

QgdsstdHEac tual Qgasstdstat ionactual Qgasstdbypa ssactual 



E) Alternatively to Calculation "C", "D" and *'F", given ongoing Qgasstdbypassactuai (from 
Calculation "H" or "I", as appropriate or "L"), TgasHEin , TgasHEout , and Tgasmixed data is available to 
the HEMD, then the actual natural gas flow rate through the heat exchanger (QgasstdHEactuai) can be 
calculated by iteration, as follows; 

^ ^Ij^ \ { \Q gasstdHEactual QgasstdbypassQctml^^'^gasmixe(\~^\Qg 

xi gasstdHEactual ~~ ^iti gasstdHEactual i^gasstdbypassactual ) 

\ TgasHEin ) 

F) Alternatively to Calculation "C", "D" and "E'\ given ongoing Qgasstdstationactuai (from 
Calculation 'T" or "S"), TgasHEin , TgasHEout » and Tgasmixed data is available to the HEMD, then the 
actual natural gas flow rate through the heat exchanger (QgasstdHEactuai) 

can be calculated by iteration, as follows; 
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G) Given ongoing QgasstdHEest (from Calculation "B"), TgasHEin , TgasHEom , and Tgasmixed data is 
available to the HEMD, then the estimated natural gas flow rate through the heat exchanger 
(Qgasstdbypassest) can be Calculated by iteration, as follows; 
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H) Given the heat exchanger bypass valve is a vertical displacement slab type (i.e. gate) valve. 
Then given ongoing APgasHE , PgasHEm ? TgasHEm , Zgas , SGgas , and BVposition data is available to the 
HEMD and converting BVposition into a heat exchanger bypass valve opening area (i.e. Asvopen) 
(Note Ybv and Cbv are based on the relationship of Aevopen to the cross-sectional area of the bypass 
piping associated with the heat exchanger bypass valve and can be obtained from lookup tables), the 
actual natural gas flow rate through the heat exchanger (Qgasstdbypassacmai) can be calculated, as 
follows; 
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I) Given the heat exchanger bypass valve is a quarter turn type (i.e. ball or plug) valve then 
APgasHE must be divided by two Then given ongoing APgasHE , PgasHEin , TgasHEin , Zgas , SGgas 5 and 
BVposition data is available to the HEMD and converting BVposition into a heat exchanger bypass valve 
opening area (i.e. ABVopen) (Note Ybv and Cav are based on the relationship of Aevopen to the cross- 
sectional area of the bypass piping associated with the heat exchanger bypass valve and can be 
obtained from lookup tables), the actual natural gas flow rate through the heat exchanger 
(Qgasstdbypassactuai) Can be Calculated, as follows; 



gasstdhypa ssactual 



0.9033 *Y,y*C,,*A,y,^^* 



gasHE 



gasHEin 

P. 



gasstd 
"TgasHEin 



0.5 
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J) Alternatively to Calculations "H" or "I", as appropriate or "K" or *'L") given on going 
Qgasstdstationactuai (from Calculation "P'* or "R") and QgasstdHEactuai (from Calculation "C" or "F"), then 
the actual natural gas flow rate through the heat exchanger (Qgasstdbypassactuai) can be calculated as 
follows; 

Qgasstdbypassactuai *~ Qgasstdstationactuai QgasstdHEactuai 



K) Alternatively to Calculation "H" or "I" as appropriate or "J" or "L", given ongoing 
QgasstdHEactuai (froHi Calculation "C" or "F"), TgasHEin , TgasHEout , and Tgasmixed data is available to the 
HEMD, then the actual natural gas flow rate through the heat exchanger (Qgasstdbypassactuai) can be 
calculated by iteration, as follows; 



Qgasstdbypassactuai 



[(« 



QgasstdHEactuai ) * 'tgasmixed \ [fij 



gasstdHEactual * 'TgasHEout I 



gasHEm 



L) Alternatively to Calculation "H" or "I" , as appropriate or "J" or "K", given ongoing 
Qgasstdstationactuai (from Calculation "P" or "R"), TgasHEin , TgasHEout . and Tgasmixed data is available to the 
HEMD, then the actual natural gas flow rate through the heat exchanger (Qgasstdbypassactuai) can be 
calculated by iteration, as follows; 

[6ga55fd!;mrio/iacA/a/ * ^gavmtrerf ] ~ [(Sgasx/f/stoftonacftw/ ~ Qgasstdbypassactuai ) 



gasstdbypassactual 



gasHEoui J 



M) Given ongoing APgasHE , PgasHEin , TgasHEin , Zgas , and SGgas data is available to the HEMD and 

that APgasHEorig 5 PgasHEinorig , TgasHEinorig , Zgasorig , SGgasorig , and QgasstdHEorig haS been previoUSly 

recorded when 100% of the natural gas flow was diverted through the heat exchanger during the 
initial calibration/programming of the HEMD. Then with the heat exchanger bypass valve 
temporarily closed to again divert 100 % of the gas flow through the heat exchanger, the periodic 
estimated gas flow rate through the compressor station (Qgasstdstationest) can be calculated as follows; 
(Note the accuracy of this calculation is dependent upon the heat exchanger fouling resistance 
remaining constant) 



Q 



gasstdstat ionorig 



gasHE 



AP 

\^ ^ gasHEorig J 



SG 



gas J 



gasHEin 



P 

\ gasHEinorig J 



(t 

gasHEinori g 

T 

\ ^ gasHEin 



( . 



' gasong 



"gas J 
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N) Alternatively to Calculation "M" or "O", given ongoing Qgasstdbypassactuai (from Calculation 
''H" or "I") and QgasstdHEest (from Calculation "B") data is available to the HEMD, then the estimated 
natural gas flow rate through the compressor station (Qgasstdstationest) can be calculated with the heat 
exchanger bypass valve in the "as found" position (i.e. without temporarily closing the bypass valve 
as in Calculation "M"), as follows; 

gasstdstat ionest " Q gasstdhypa ssactuai QgasstdHEes t 



O) Alternatively to Calculation "M" or "N", given ongoing QgasstdHEest (from Calculation "B"), 
TgasHEin , TgasHEout , ^ud Tgasmixed data is available to the HEMD, then the actual natural gas flow rate 
through the heat exchanger (Qgasstdstationest) can be calculated by iteration, as follows; 



C gasstdstationest 



=0 



gasstdHEest 



*r 

gasstdstationest * gasmixed 



gasstdHEest ^ gasHEout 



^gflsHEin 



P) Given on going APgasmeteKstation) , Pgasmeter(station) , Tgasmixed , Zgas , SGgas , and Station gaS flOW 

orifice meter opening area (i.e. Aonficestn) , data is available to the HEMD (Note Yonfice and Confice are 
based on the relationship of Aorificesm to the cross-sectional area of the piping associated with the 
orifice meter and can be obtained from lookup tables). Then the actual natural gas flow rate through 
the compressor station (Qgasstdstationactuai) can be calculated as follows; 



e= 1 2775 * F * C ^ A * 
gasstdstat ionactual ^ orificestn orificestn orificestn 



(station ) 



gas J 



\ 



gasmeter (station ) 
P 

y * gasstd J 



0.5 



' gasstd 



Q) Altematively to Calculation "P", "R" or "S"„ given on going Qgasstdbypassactuai (from 
Calculation "H" or "I", as appropriate or "K") and QgasstdHEactuai (from Calculation "C" or "E") data 
is available to the HEMD, then the actual natural gas flow rate through the compressor station 
(Qgasstdstationactuai) Can be Calculated as follows; 

Q gasstdstat iomctml Q gasstdhypa ssactuai QgasstdHEac tual 



R) Altematively to Calculation "P", "Q" or "S", given ongoing QgasstdHEactuai (from Calculation 
"C" or "E"), TgasHEin , TgasHEout , and Tgasmixed data is available to the HEMD, then the actual natural 
gas flow rate through the heat exchanger (Qgasstdstationactuai) can be calculated by iteration, as follows; 
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' gasstdstationactual iC-gasstdHEactual 



QgasstdHEactual 



gasstdstationactual 



S) Alternatively to Calculation "P", "Q" or "R", given ongoing Qgasstdbypassactuai (from 
Calculation "H" or *T', as appropriate or "K"), TgasHEm , TgasHEout , and Tgasmixed data is available to 
the HEMD, then the actual natural gas flow rate through the heat exchanger (Qgasstdstationactuai) can be 
calculated by iteration, as follows; 

[^g<U5ftbi:a/ionacA«i/ * ^osmoccrrf ]~" [(^gosslb^to/io/iac/iia/ ~" Qgasstdbypassactuai )* '^gasHEota ] 



gasstdstationactual 



-a 



gasstdbypassactual 



gasHE'm 



T) Given on going PgasHEin , TgasHEm , Zgas , SGgas , and QgasstdHEactual (from Calculation "C", "D'\ 
"E" or "F") data is available to the HEMD and that APgasHEong , PgasHEinorig , TgasHEinong » Zgasong , 
SGgasorig » and QgasstdHEorig has been previously recorded when 100% of the natural gas flow was 
diverted through the heat exchanger during the initial calibration/programming of the HEMD, then 
the expected pressure differential across the heat exchanger (APgasHEexpect) can be calculated as 
follows; 



A P = AP * 

g(isHE exp ect gasHEorig 



{Q. 



gasstdHEac tual 



SG 



gasHEinori g 
P 

\ ^ gflsHEin J 



gasHEin 

T 

\^ gasHEinori g J 



\ QgasstdHEorig J 

U) Given on going APgasHE and APgasHEexpect (from Calculation 'T") data is available to the 
HEMD, then the change in pressure differential across the heat exchanger (APgasHEchange) can be 
calculated, as follows; 



AP = AP - AP 

^gasHEchang e ^gasHE ^ gasHE cxp ect 



If APgasHEchange > 0 (zero) then an increase in pressure differential across the heat exchanger has 
occurred. Conversely, if APgasHEchange < 0 (zero) then a decrease in pressure differential has occurred. 
In any case the HEMD will then send either an open or close (BVopen or BVdose) command signal to 
the heat exchanger bypass valve operator to adjust the bypass valve position (BVposition) and thereby 
adjust the natural gas flow through the exchanger in order to accommodate the new conditions and 
obtain the new optimum balance between gas cooling and heat exchanger pressure loss. The HEMD 
will record any increase in pressure differential across the heat exchanger (APgasHEchange) and flag the 
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incident for possible fixture maintenance action as the increase may be due to internal fouling of the 
exchanger tubes. 



V) Given TgasHEinong , TgasHEoutorig , wfcgasHEorig 5 Cpgasorig 5 AhEcxI , Forig 5 and LMTDorig, data is 

available to the HEMD after initial calibration/programming of the HEMD, then the actual original 
overall heat trmisfer coefficient (UHEdirtyorig) with fouling resistance (Rfouiing) (i.e. dirty) included can 
be calculated, as follows; 



U 



HEtlirtyorig 



^^^asHEorig * Cpgasorig 

*LMTD^^. 



* gasHEinorig ^ gasHEotitortg 



) 



W) Given on going TgasHEin , TgasHEout , ^&gasHE , Cpgas , AHEext , F, and LMTD, data is available to 
the HEMD, then the actual overall heat transfer coefficient (UnEdirtyactuai) with fouling resistance 
(Rfouimg) (i e. dirty) included can be calculated, as follows; 



U 



HEdirtyactual 



\asHE * ^pgas 



K^HEe«*P*LMTD 



[t 



* ' ^ gasHEin 'TgasHEout 



] 



X) Given dtube , Dtube , ^tuhe 9 hgasorig and hairorig data is available to the HEMD after initial 
calibration/programming of the HEMD, then the original overall heat transfer coefficient 

(UHEcieanorig) with no fouliug resistance (Rfouiing) (i.e. clean) can be calculated, as follows; 

-1 



U 



HEcleanorig 



1000* 



1 



tube 



*h 



gasorig 



1 



2 ^ '^tube 



tube 



''tube 



1 



^ 



\^airorig J 



Y) Given on going dtube > Dmbe j ktube > hgas and hair data is available to the HEMD, then the actual 
overall heat transfer coefficient (UnEcieanactuai) with no fouling resist^ice (Rfouimg) (i.e. clean) can be 
calculated, as follows; 

7 \ ( ^ 



u 



HEcleanactual 



1000* 



gos 
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2 * ^"^^ 



,^^^£tube_ 



*'tube 



tube 
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Z, Given U„.,^. ^ Uh.^-. calculate. "V and "X" J^" 
L^^r^ flow rates trough me hea. exchanger, da. . avaiUM. .0 *e *e 
ori^na. fouling r^istance (R..^^ ^ occurred in Are hea. exchanger can be calcula«d as 



follows; 

1 1 



R 



■foulingorig Tf IT.,.. , 



J T T rfi-ATTi calculations "W" and "Y" based upon the 

AA) Given on going UHEdirtya-^tuai and UHEcieanactua. (from calculations vv ^^^^^^ 

Jenatura. gas nrassflowrarcsd^ughdrehear exchanger, da.isavai,ah.e.*e^^ 

■ . (V ,^ that has occurred in the heat exchanger can be calculated 

the actual fouUng resistance (Rfouimgactuai) tnat nas oo^micu 



as follows; 

1 L__ 



BB, Given on going R.,^. and Rr-„^, (fix,m calculations and "AA", dam is available 
.0 fte HEMD, to Ure actual change in fouling resistance (R«^ *at has occurred rn the heat 

exchanger can be calculated, as follows; 

Rfoulingchcnge = ^foulingi^tml ' ^foulmg»rig 

CO Given on going R....^. (from calculation "BB") and AP,« (from calculation "S") 
data is available to the HEMD, then it is possible to detect whether or not the difference in fouling 
resistance (Rfouiing) is internal or external, as follows; 

. n r^ero^ then fouUng has been removed from the heat exchanger 
If in the first case Rfouiingchange < 0 (zero), tnen louung n<«, . ^ . ,• 

tubes If in the second case R-„.h„. - 0 (zero, then no change has occurred m dre fouhng 
resistanceofdrehea. exchanger n,bes. If in the fl,irdcaseR«.«>0(zero,andAP^B.^>^ 

(zero, then an increase in internal fouling resistance of the heat exchanger tubes has occurr^. If m 
the .burdr case R«.^ > 0 (zero) and AP.„^ - 0 (zero, then an increase in extenra. 
fouling resistance of dre hea. exchanger tabes has occurred. In the first case the abthty o he h^ 
exchanger .o cool the natural gas has been enhanced. In the second and dW cases dre abth^- of dre 
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heat exchanger to cool the gas has been diminished. The HEMD will calculate and quantify the 
impact of the change in fouling resistance on the heat exchanger outlet gas temperature using some 
or all of the equations found in Paragraphs 17 to 32. The HEMD, will then send either an open or 
close (BVopen or BVciose) command signal to the heat exchanger bypass valve operator to adjust the 
bypass valve position (BVposition) and thereby adjust the natural gas flow through the exchanger, in 
order to accommodate the new conditions and obtain the new optimum balance between gas cooling 
and heat exchanger pressure loss. The HEMD will record any increase in fouling resistance of the 
heat exchanger tubes (ARfouUngchange) and flag the incident for possible future maintenance action. 

DD) Given on going TgasHEin and TgasHEout data is available to the HEMD, then the temperature 
differential across the heat exchanger (ATgasHE) can be calculated as follows; 

AT =T 

gasHE ^ gasHEin ^ gasHEout 

EE) Given on going TgasHEin and Tairamb(HEin) data is available to the HEMD, then the temperature 
differential between the hot heat exchanger inlet natural gas temperature and cool exchanger inlet 
air temperature (ATgas/air) can be calculated as follows; 

gas I air ~ TgasHEin ~' airamhiHEin) 

FF) Given on going APgasHE , ATgas/air , and Qgasstdstation data is available to the HEMD, the 
optimum heat exchanger pressure differential (APgasHE{optimum)) and thereby bypass valve opening 
area (ABVopen(optimum)) is determined by the control algorithm. The difference (AevopenCchange)) between 
the optimum opening area and the current bypass valve opening area (ABVopen(cuiTcnt)) can be 
calculated, as follows; 
A -A -A 

BVopenichange) "^BVopenioptimum) ^BVopen{cun&it) 

If ABVopen(change) > 0 (zcro) then the HEMD will convert ABVopen(change) into the corresponding 
optimimi heat exchanger bypass valve position (BVposition(optimum))- The HEMD will then send an 
open command signal (BVopen) to the heat exchanger bypass valve operator to open the bypass 
valve until the optimum opening position has been achieved. If ABVopen(changc) < 0 then the HEMD 
will again convert ABVopen(change) into the corresponding optimum heat exchanger bypass valve 
position (BVposition(optimum))- The HEMD will then send a close command signal (BVdose) to the heat 
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exchanger bypass valve operator to close the bypass valve until the optimum opening position has 
been achieved. 



GG) Given Qgasstdstationactuai (from outside sources), Qgasstdstationest (from Calculation "0") and 
QgasstdHEnest (from Calculation "B"), then the actual natural gas flow rate through the heat exchanger 
(QgasstdHEactuai) can be calculated, as follows; 



Q gasstdHEactuai QgasstdHEest * 



^ gasstdstationactual 
\ Qgasstdstationest ^ 



11 The controller 46 is preferably operated to achieve optimum flow energy savings. This may 
be achieved by reference to the graphs of Figs. 14 to 17, or their mathematical equivalent stored in a 
computer. 

78 In Figs. 14 and 15, the curves are based on a 30 MW compressor station with heat 
exchanger. The lower curve on this graph represents the incremental compressor power required to 
overcome the pressure loss associated with flowing gas through a heat exchanger. The far left of 
this curve indicates the maximum compressor power required to overcome the pressure loss 
associated with flowing all of the natural gas from the compressor through the heat exchanger. The 
right end of the curve shows that no compression power is required when there is no flow through 
the exchanger 



79 The upper curve depicts the power saving that results from transporting cooled natural gas 
tluough a pipeline, when all exchanger fans are operating, based on the cooling available for a 
typical Spring or Fall day. The upper curve drops off slowly because the amount of heat transferred 
(from the natural gas to the air in this case) is not directly proportional to the gas flow rate through 
the exchanger. For example when gas flow rate through the exchanger is reduced to one half the full 
gas flow rate, the amount of heat transfer equals approxunately 78 % of the fiiU flow heat transfer 
(for an average sized mainline transmission heat exchanger) (see equations in Paragraphs 28 and 
29). In this case although less natural gas is being cooled at the reduced flow rate, the gas is being 
cooled to a greater degree (i.e. closer to the ambient air temperature). The net effect is a non- 
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proportional change in the amount of heat transferred with a reduction in gas flow rate through the 
exchanger, as depicted in this graph, 

80 Upon comparison of the two curves it can be seen that the slope of the upper curve initially 
changes at a slower rate than that of the lower curve. It is this phenomena that provides the basis for 
optimizing exchanger performance. The difference between the two curves represents flie net power 
saving from cooling the natural gas. For example 80 % of the gas flow rate through the exchanger 
requires an extra 1400 kW of compressor power (to overcome exchanger pressure loss), but results 
in a (gross) savings of 2800 kW due to a lower pressure loss along the downstream pipeline 
(resulting from lower temperature inlet gas). This results in a net power saving of approximately 
1400 kW. 

81 As shown in Fig. 15, by comparing the two curves depicted in the graph of Fig. 14, it can be 
seen that the greatest net power saving occurs when approximately 50 % of the natural gas flows 
through the exchanger (i.e. optimum gas flow rate occurs when approximately 50 % of the gas flow 
rate bypasses the exchanger) and results in a net power saving of approximately 1750 kW. 

82 Referring to Fig. 16, which is a graph showing incremental power saving curves based on 
the cooling available for an average winter. Spring or Fall, and Summer day and the incremental 
compressor power curve, the lower power saving curve represents the impact of the heat exchanger 
based on a Summer day. The greatest difference between this curve and the incremental compressor 
power curve occurs when approximately 34 % of the natural gas flows through the exchanger (i.e. 
optimum gas flow rate occurs when approximately 66 % of the gas flow rate bypasses the 
exchanger) and results in a net power saving of approximately 600 kW. The center curve is the 
same one depicted in Fig. 14 and is based on the cooling available for an average Spring or Fall day. 
As stated previously the greatest net power saving occurs when approximately 50 % of the gas 
flows through the exchanger. 

83 The upper curve represents the impact of the heat exchanger based on an average winter 
day. The greatest difference between this curve and the incremental compressor power curve occurs 
when approximately 60 % of the natural gas flows through the exchanger (i.e. optimum gas flow 
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rate occurs when approximately 40 % of the gas flow rate bypasses the exchanger) and results in a 
net power saving of approximately 1950 kW. 

84 Referring to Fig. 17, it can be seen that a relationship exists between the temperature 
differential between the heat exchanger inlet gas and ambient air (i.e. ATI, AT2 and ATS) and the 
optimum gas flow rate through the exchanger. This in turn can be translated into a relationship 
between the temperature differential between the exchanger inlet gas and ambient air and the 
optimum heat exchanger gas pressure differential which is the basis for the heat exchanger control 
algorithm. This relationship is unique to every heat exchanger installation and is based upon 
numerous parameters including heat exchanger design (i.e. total cooling surface area, and number, 
configuration, length and diameter of the cooling tubes, and overall heat transfer coefficient), heat 
exchanger inlet gas temperature and pressure, ambient air temperature, the specific heat capacity of 
the natural gas and air, total gas flow rate fi:om upstream compressor(s), maximum air flow rate 
through the exchanger, compressor performance and efficiency, downstream pipeline 
characteristics (i.e. pipe size, wall thickness, burial depth, etc.) and soil conductivity. 

85 The net power saving resulting fi-om optimizing natural gas flow rate through the exchanger 
results in a reduction in the power to gas throughput ratio and can improve the efficiency of a 
pipeline in one of two ways. The net power saving can either be utilized to move the same amount 
of gas through a pipeline with less compressor power or to move more gas through the pipeline with 
the same amount of power as is currently required. The reduction in the power to gas throughput 
ratio also results in a significant reduction in CO2 emissions per given gas throughput. 

86 The following nomenclature is used in the equations presented in this disclosure: 



Abundie - heat exchanger bundle area available for air flow (i.e. total bundle area less area taken up by 
tubes), 

Aevopen = beat exchanger bypass valve opening area based upon BVposition» ni^ 

Ahecxi = total bare tube external surface area of heat exchanger, 

AorificeHE = hcat cxchanger gas flow meter orifice opening area, 

Aorificestn ~ statlon gas flow meter orifice opening area, 

Atube = cross-sectional area of heat exchanger cooling tube, 

BVciose - heat exchanger bypass valve close command 

BVopen = heat exchanger bypass valve open conmiand 

BVposirion = heat exchanger bypass valve actual position, (0 - 90 degrees for quarter turn valves and 0 - 

100 percent of total travel for linear actuated valves) 
Cpair = specific heat capacity of air at constant pressure, kJ/kg-^K 
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Heal CACIlallgCI L'UUllllg lUUC lllolUC UlalllClCly 111 


J-'tube 


hpat pxrhariper ronlitip tiihe oiitJ?ide diameter m 








D'Arcy friction factor of heat exchanger cooling tubes 


F = 


correction factor when using other than a counterflow double pipe type heat exchanger (i.e. 
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heat transfer coefficient of natural gas, W/m '^K 
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kilnPascak T 1 000 Newton/m^^ 


kW = 


kilowatt (kilojoule/s) 
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adiabatic exnonent of natural &as c«»ac /c«««e 




thermal conductivity of air, W/m °K 


^as ~ 


thermal conductivity of natural gas, W/m °K 


ktube ~ 


thermal conductivity of heat exchanger cooling tube, W/m-^K 


Lpipe ~ 


length of pipeline, km 


Ltube ~ 


length of heat exchanger cooling tubes, m 
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heat exchanger cooling air mass flow rate, kg/s 




natural gas mass flow rate, kg/s 


'^gasHE = 


natural gas mass flow rate through heat exchanger, kg/s 




natural gas mass flow rate through compressor, kg/s 


'"^gaspipe 


natural ffa^ mass flow rate through niTieline kff/s 
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adiabatic con^ressor efficiency 
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ambient air pressure, kPa(a) 
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T>ir»pliTip miflpf natural crac nrpssiirp VPa/^a^ 


P d ~ 


natural gas pressure at standard conditions, 101.325 kPa(a) 


Pgascompdis ~" 


compressor discharge natural gas pressure, kPa(a) 


P = 

gascompsuc 


con^ressor suction natural gas pressure (sPgas compdis -APgasHE), kPa(a) 


PgasHEin " 


heat exchanger inlet natural gas pressure, kPa(a) 


P gasineter<station) ~ 


natural gas pressure at compressor station gas flow meter, kPa(a) 


APgasHE ~ 


natural gas pressure differential across heat exchanger (i.e. pressure loss), kPa 


^Pgasmeter(HE) ^ 


natural gas pressure differential across heat exchanger gas flow meter, kPa 


^Pgasmeter(station) ~~ 


natural gas pressure differential across compressor station gas flow meter, kPa 
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heat (ertfiTPv) transfer hv heat exrhanper IcT 


Qgasstdbypass 


calculated flow rate of natural gas through heat exchanger bypass valve at MSC, m^/s 


(JnaectriHF ~ 


calculated flow rate of natural gas through heat exchanger at MSC, m^/s 


Qgasstdstation 


calculated flow rate of natural gas through compressor station at MSC, m^/s. 


R 


universal gas constant equals 8.314 kJ/kmol-^K 


■^fouling 

S = 


heat exchanger tube fouling resistance factor , m^-°K/kW 
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SGgas ^ 


specific gravity of natural gas 


'^airainb(HEin) ~ 


ambient air temperature (i.e. heat exchanger air inlet ten:q)erature), 
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heat exchanger outlet air ten^erature , ®K 
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compressor discharge natural gas temperature, **K 
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■*• gESCompsuc 


compressor suction natural gas temperature, °K 
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XgasHEout ~ 


heat exchanger outlet natural gas temperature, **K 




pipeline inlet natural gas temperature , °K 


Tgasmixed 


natural gas temperature resulting from the mixture of the heat exchanger outlet gas (TgasHEout) 
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pipeline outlet natural gas temperature , °K 


gasstd 


natural gas temperature at standard conditions, 288.15 °K (15 °C) 
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natural gas temperature differential between the heat exchanger inlet gas and ambient air, °K 
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^ gascomp 


natural pas temtierature differential between the comnressor inlet and outlet cas °1C 


ATgasHE ^ 


natural ffas tpmneratiire diffprpntial betwppn the hpat exchanger inlpt and outlet pas ®1C 




UVClall iiCal liaiioiCl V,^UCiiivivill Ui liCal CAl/iiailgvi, ivvv/Xii iV 




uynauiic viscosiiy oi au, Kg/m'S (^inoic i Kg/m'S — iukju ^cmipoisc^ 




dynamic viscosity of natural gas, kg/m-s (Note 1 kg/m-s = 1000 Centipoise) 


■117 

w - 


Watt 




ilCdl CA^UaiigCl Uypabo VaiVC IlCl CApallblUIl iav/lUI ^iiUili voUICd^ 


Yorifice ^ 


gas flow orifice meter net expansion factor (from tables) 


Zgas ~ 


natural gas compressibility (calculated based upon gas components, pressure and 




temperature) 




compressor discharge natural gas con^ressibility (calculated) 


Zgascompsuc 


= compressor suction natural gas con:q)ressibility (calculated) 


Zgasinlet ~ 


pipeline inlet natural gas conpressibility (calculated) 


Zgasoutlet ~ 


pipeline outlet natural gas con^>ressibility (calculated) 
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